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ABSTRACT 

In this study, the offset strip fin array in a liquid flow through module was 
examined. The study used a model of the heat exchanger which is ten times larger than 
the actual model used in electronics cooling. Water and Brayco Micronic 889 
Polyalfrfiaolefin were used as the working fluids in order to determine the Reynolds 
number and Prandtl number effects on the Colburn j factor. Also, flow visualization 
experiments were performed to determine the flow characteristics of the offset strip fins 


for lamiruir flow. 
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L INTRODUCTION 


A. ELECTRONIC COOLING 

The trend in electronic packaging has been towards building more electronic 
circuits into one very small chip. Since the invention of the silicon integrated circuit, 
there has been an increase in the number of components per chip. This increase in circuit 
integration continues to be driven by the need for increased processing speed. With each 
increase in the scale of integration of microelectronic circuits more effective cooling 
systems are required. 

Currently, chip heat fluxes are approaching 50W/cm^ and heat fluxes 
of 200 W/cm^ are projected [Ref 1]. Also, it is known that chip temperatures must be 
maintaine d within prescribed operating limits below 8S° C. Due to the large heat fluxes 
and small areas associated with current aisd projected very large scale integrated circuits 
and the maximum allowable junction temperature requirement, it could be necessary to 
cool the system by liquid cooling and to enharKe heat transfer by using extended 
surfaces. 

B. UQUID COOLING METHOD 

At lower air velocities, typical heat flux limits for packages cooled by 
conventional forced convection are less than 1 W/cm^ [Ref 2]. Because of the hi^ heat 
fluxes of today's technology, it is evident that air cooling could not provide an optimal 
solution to the cooling problem. As a result of the higher density and thermal 
conductivity of liquids in comparison to air, it is possible to enhance heat transfer 
coefficients by using liquid cooling. 

1. Direct Liquid Cooling 

Direct liquid cooling maintains physical contact between the coolant and the 
electronic components. The coolant should have dielectric characteristics and good 
chemical compatibility so as not to effect the circuit performance. 



2. Indirect Liquid Cooling 

Indirect liquid cooling cools the electronic components with attachment of a 
liquid cooled plate/heat exchanger to the electronic package (See Figure 1). This cooling 
technique requires the addition of a liquid circulating piunp. The significant advantage of 
this cooling method is the reduction in the thermal resistance to heat transfer. 

Modules are usually enclosed in protective covers to which fins may be attached 
to enhance heat transfer to an external coolant. 



C. OFFSET STRIP FINS 

1. Characteristics of Offeet Strip Fins 

Offset strip fin is one of the most widely used fin geometries in compact heat 
exchangers. The cross section of the fin is rectangular. The fin is cut into small strips of 
length I in the flow direction. The other variables are the fin spacing s, fin height h, and 
fin thickness t. The offset is about 50 % of the fin pitch (See Figure 2). 






The offset strip fins have high heat transfer coefficients. This is mainly caused by 
the developing laminar boundary layers. The interrupted geometry of the fins prevents 
the velocity and temperature boundary layers from getting excessively thick. Since the 
form drag increases due to the finite thickness of the fins the firiction factors are also 
high. In this study the offset strip fin array in a liquid flow through module (LFM) used in 
electronics cooling is examined. The study uses a model of the heat exchanger which is 
ten times larger than the actual model. A typical liquid flow through module is shown in 
Figure 3. 



2. Analytical and Numerical Correlations for Offset Strip Fins 

Weiting [Ref 3] correlated available experimental heat transfer and flow fiiction 
data for 22 offset strip fin surfaces using air as the working fluid. Weiting defines the 
Reynolds number based on the hydraulic diameter as follows; 
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. Liquid Flow Through Module. 







Re=^ 


Equation 1.1 


Here U is the mean velocity of the air in the channel, v is the kinematic viscosity 
of the air and D* is the hydraulic diameter (Di, = 4AclP) where Ac represents the 
cross-sectional area and P represents the perimeter of the channel. Weiting’s correlations 
were as follows; 


For Reo i 1000, 


/ = 0.483(//Da)-^ >“(a*)-o 


For Reo ^ 2000, 

f= 1.136(//Z)*)-«’**(8/D*)» 

/ = 0.242(//D*)^ 322(5/£,^)0.089^g^.368 


Equation 1.2 
Equation 1.3 


Equation 1.4 
Equation l.S 


Here / is the strip length, S is the fin thickness, D* is the hydraulic diameter of the 
passages and a* is the ratio of the width to height of the passage. For Reynolds munbers 
in the laminar range he stated that the Colburn j factor was only a function of Reynolds 
number and //Da or a* . He considered the correlation to be accurate within ± 10 % , 
though some points deviated by forty percent. 

Manglik and Bergles [Ref 4], gathered all the previously taken data for the 
offset strip fins and developed a correlation for the Colburn j factor which can be used 
for both the laminar and turbulent regions. Their definition of Reynolds number and 
hydraulic diameter was the same with Weiting’s definition. 




Equation 1.6 



In their correlation, s is the fin spacing, h is the fin height, t is the fin thickness 
and 1 is the fin length. They state that this correlation is applicable for Prandtl numbers in 
the range of 0.5 to 15. They estimated their correlation to be accurate within ± 20 %. 

The correlations for the transition tfom laminar to turbulent region were 
developed by H. M. Joshi and R. Webb [Ref. 5]. They used a numerical solution for the 
Nusselt number and the fiiction factor for a laminar flow and a semi-emprical solution 
for the turbulent region. 

For Rez) < 1000, 

/ = 0.53{//Z)*)-<'‘5(^/;,)-oo2/jg^.5o Equation 1.7 

For /?ez)^/?ez)+ 1000, 

/ = 0.21(//DA)-^24(,/£,^)OO2;jg-o4o Equation 1.8 

Several attempts have been made to predict analytically the heat transfer and 
flow fiiction characteristics of offset strip fins. One of the first attempts was made by 
Kays [Ref 6]. He used Pohlhausen l amin ar boundary layer solution for a flat plate to 
calculate the heat transfer from an offset strip fin. For fiiction, he used the Blasius 
laminar boundary layer solution for a flat plate. His Reynolds number was not based on 
the hydraulic diameter but on the length of the flat plate {Rei = U ■ l/v). His results were 
as follows; 

i = 0.6MRe^^ Equation 1.9 

f=Y+^ .328Re^^ Equation 1.10 

Here 1 is the fin length, 5 is the fin thickness and Cq is the drag coefficient. 
Kays suggested Cd= 0.88 based on the potential flow normal to a flat plate. In this model 
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it is idealized that the velocity and temperature boundary layers are destroyed in the wake 
region. Although this model is at best a first order approximation, it does not provide the 
detailed insight into the systems . Patankar and Prakash [Ref 7], predicted the j and f 
factors for the offset strip fins numerically by a finite difference method. They 
considered the laminar boundary layers on each strip fin They first assumed that strip 
fins have infinite height and zero thickness. Then they extended the analysis for finite fin 
thickness. They observed that at low Reynolds Numbers or low 6/s, the recirculation 
zone behind the trailing edge was small. At high Reynolds numbers or high 8/s, this 
recirculation zone extended from the trailing edge to the next strip. (See Figure 4). 

Although a large amount of experimentation and numerical predictions have 
been made, further studies and improvements are needed to accurately understand the 
performance of the offset-strip fins. 




Figure 4. Flow Patterns Over Offset Strip Fins At Different 8/s for Rez)=2000, Based 
On A Numerical Solution, From Ref 7. 




D. 


OBJECTIVES OF PRESENT STUDY 


The first objective of this study was to investigate the heat transfer 
characteristics of a liquid flow through module which is actually a ten times scaled up 
model of the fluid circulation passages of an SEM-E (Liquid Flow Through Module). By 
using water and Brayco Micronic 889 (Castrol) Polyalphaolefin (PAO) as the working 
fluid, experiments were performed to determine the Reynolds number and Prandtl 
number effects on the dimensionless heat transfer coefficients. 

The second objective was to visualize the flow over the offset strip fins which 
were used to enhance the heat transfer. Dye injection was applied to observe the variation 
in the velocity boundary layer across the fins and the recirculation of the flow behind the 
trailing e<^e of the fins for varying flow rate of the coolant. 



IL EXPERIMENTAL APPARATUS 


The experimental apparatus included a system to circulate the coolant through 
the test section and another system to collect the necessary data to be able to calculate the 
heat transfer and flow characteristics of the test section. 


A. COOLANT CIRCULATING SYSTEM 

A schematic of the coolant circulation system is shown in Figure 5. 



1. Liquid Tank 

In this study first water and then Brayco Micronic 889 (Castrol) Polyalphaolefin 
(PAO) were used as the working fluids. In both case, the coolant was provided by an 
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Endocal RTE-5 circulating bath. The Endocal RTE-5 has a heater/refngeration unit 
which enables the regulation of the temperature of the liquid in the bath between -30 to 
100 °C. Eturing the experiments this temperature was maintained constant at around 
20“C. Since the volume of the circulating bath was small, it was noted that as the 
Reynolds number exceeded about 650 for water and 60 for Brayco Micronic 889 
Polyalphaolefin, the temperature of the liquid in the bath was effected by the hotter liquid 
at the exit of the heat exchanger which was circulated back to the tank. 

2. Pump 

A Cole Parmer positive displacement gear type pump was used to circulate the 
coolant through the heat exchanger. The positive pressure fluid surge volume to the 
pump was provided by the heater/refrigeration unit of the liquid bath. 

The pump was driven by a DC motor. The speed control unit of the DC motor 
provided varying flow rates and Reynolds numbers for varying pump speeds. For the heat 
transfer experiments the pump speed was varied between 10-70 % of its maximum value 
which corresponded to the Reynolds numbers in the range of 50-650 for water and the 
Reynolds numbers in the range of 10-80 for Brayco Micronic 889 Polyalphaolefin. 

It was noted that for small pump speeds the time required by the test section to 
reach its steady state condition was around 60 minutes. The time required to reach the 
steady state condition was decreasing with increasing pump speed (increasing Reynolds 
number). 


3. Flowmeter 

Following the pump, coolant was passed through an Omega FTB-102 turbine 
flowmeter to measure the flow rate. As the flow rate varied the voltage across the AC to 
DC signal converter which was connected to the flowmeter also changed. These voltage 
outputs were than used to calibrate the flowmeter. 
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The flowmeter was calibrated for both of the coolants. The calibration of the 
flowmeter for water was made in two different ways. The first one consisted of a 
comparison of the average voltage output to the average time to collect a specific volume 
of water. The second method was based on the collected mass of the water in a specific 
time for varyii^ pump speeds. Even though the second method produced a more linear 
calibration curve than the first one, the flow rate measurements, especially in the range of 
middle Reynolds numbers (200-600), were very close to each other. The calibration of 
the flowmeter for the PAO was based on the first method. All flowmeter calibration 
curves were presented as Appendix A. 

4. Heat Exchanger 

The flowmeter outlet was connected to a manifold and divided into three flow 
paths before entering the offset strip fin heat exchanger. 

a. Offset Strip Fins 

The test section wliich includes the base plate and the offset strip fin 
assembly was first constructed by Lt Jeffrey Masterson [Ref 8]. A picture of the test 
section is shown in Figure 6. The fin assembly was a 10 times scaled up model of the 
horizontal passages of a liquid flow through module used in electronics cooling. The base 
plate consisted of three parts. (The inlet plenum, outlet plenum and the middle part with 
fins). The inlet and outlet plenum were made from two pieces of a 12.7 mm aluminum 
6061 alloy plate. They were milled down to the base plate thickness of 10.16 mm. The 
middle part of the base plate, which included the offset strip fins, was also made from 
aluminum 6061 alloy. Instead of welding the fins onto the base plate, the base and the 
fins were milled from a solid piece, not to leave weld deposits and affect the fin 
geometry. A 2.54 cm thick aluminum plate was cut into 31.75 mm wide strips. These 
strips were then milled to the width of the plate. The fins were cut 15.24 mm deep into 
the aluminum and 15 rows were obtained. Each row had 19 fins where the spacing 
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Figure 6. Test Section. 


12 







between the fins was 11.76 mm. A 6.64 nun offset was allowed on one end that when 
alternated created the offset pattern of succeeding fin rows. On each end of the base plate 
an additional 12.7 mm offset of base plate thickness was milled to attach the plexiglas 
cover [Ref 8]. A schematic of the finned section is shown in Figure 7. 



The inlet plenum, firmed section and the outlet plenum were coimected to each 
other by means of two threaded rods. The holes to accommodate the threaded rods were 
drilled in the lengthwise direction through each fin row and plenum piece. Rubber gasket 
was placed between the first row and the inlet plenum sections and between the last row 
and the outlet plenum sections to reduce the longitudinal conduction losses fi'om the 
firmed section to the plenum The final assembly was held and tightened by using nuts at 
each end of the rods [Ref 8]. 
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figure 8. Overall Dimenaons of The Test Section. 


Since the water left some deposits on the base plate and offset strip fins, 
the test section was cleaned and sandblasted and prepared for the present study. 

b. Plexiglas Cover 

The construction of the heat exchanger was completed by enclosii^ the 
base plate and fin assembly by a plexiglas cover. Since the original plexiglas cover was 
deformed a new one was fabricated. The top was cut to the dimensions of the base plate 
(654x260.5 mm) and the sides were cut to the height of the fins (15.24 mm). 

Before attaching the plexiglas cover to the base, all the necessaiy holes 
were drilled into the plexiglas cover. Three inlets and three outlets were drilled into the 
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axial ends of the plexiglas cover. Inlets and outlets consisted of 3.175 nun pipe threads 
drilled and tapped into the plexiglas wall with 6.3S nun tygon adaptors on the opposite 
end of the wall. One inlet and one outlet were aligned on the longitudinal axis of the test 
plate. The two other inlets and outlets were placed 37.9 nun from the outer edges of the 
plate. Valves were placed on the two other outlets so they could be closed when 
operating the test section as needed or opened to drain the system after each experiment. 
During the first runs these outlets with valves were kept closed and the inlet to the 
circulating bath was provided by only one outlet of the heat exchanger. But since the 
pressure in the heat exchanger was high especially for Re > 500, one of the outlet valves 
was also kept open v«Mch provided the inlet to the circulating bath from two flow paths. 
This decreased the high back pressure in the heat exchanger. 

Another step before attaching the cover to the base plate was to drill the 
dye injection and temperature probe points into the plexiglas cover. The details about the 
dye injection points will be given in the next section. 

Four temperature probe holes were placed above both the inlet and outlet 
plenum regions to measure the inlet and outlet temperatures of the coolant. These holes 
were very tiny and filled with silicone before feeding the thermocouples through the 
holes. 

Prior to attaching the cover to the base a flow straightener was also fixed 
between the inlet temperature probe and dye injection holes. The original flow 
straightener which was devised using 38.1 mm long straws and a plastic mesh was glued 
to the plexiglas cover. The purpose was to provide a uniform flow velocity profile before 
the first fin row. The distance between the first fin row and the edge of the flow 
straightener was about 28 mm. After drilling all the necessary holes and fixing the flow 
straightener, the sides were attached to the top piece using acrylic cement and secured 
with screws. During the first runs leakage was noted from the temperature probe holes 
and the side walls of the plexiglas cover. Silicone was run along the sides and refilled 
into the temperature probe boles. When the heat transfer experiments were completed a 
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crack was noted on the top piece of the plexiglas cover. Acrylic cement was syringed into 
the crack to fill the gaps. 

The heat exchanger was the last element of the coolant circulating system. The 
outlet of the heat exchanger was connected to the refrigerated bath and the circulation of 
the coolant vras completed. Flexible Tygon tubing was used to connect the inlet and 
outlet parts of the coolant circulating system. Clamps were affixed to tighten the 
connections. 

B. DYE INJECTION SYSTEM 

A schematic of the dye injection system is shown in Figure 9. 



Figure 9. Dye Injection System. 
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Five dye injection points were placed in the plexiglas cover. They were 1.6 mm 
holes that penetrated the plexiglas cover just before the first row fins and after flow 
straightener assembly. Five adaptors were placed into these holes. Inside the plexiglas 
cover, stainless steel surgical tubing was passed through the adaptors. This tubing was 
bent into a 90 degree bend and the ends were smoothed and shaped to a cone in order to 
direct the dye flow between the fins. The total length of a steel tubing was about 50 mm. 
The other ends of the adaptors which were left on the outside of the plexiglas cover were 
coimected to a manifold by using 190 mm long plastic tubes. The manifold was used to 
distribute the dye in the reservoir to the desired injection tubings. The reservoir was a 
small plastic bottle and was held about I meter above the test section to (xrovide a 
positive pressure. A small syringe also served as an ink reservoir. 

C. DATA COLLECTING SYSTEM 

A heater was used to heat the test section and the temperatures were measured 
by means of thermocouples. Thermocouples were connected to a data acquisition system 
which also measured the required voltages. Flow visualization data was recorded by 
using a 35 mm camera. 

1. Heater 

Printed circuit boards usually dissipate constant heat flux. This was simulated 
by using an electrical heater ^ch was a Minco foil backed patch. The heater pad was 
25.4 cm by 30.48 cm, resulting in a heater surface area of 774.2 cm^. 

The heater pad was fixed to the bottom of the aluminum plate with pressure 
sensitive adhesive. The pad covered the finned area section except 3.24 mm on the sides 
and 9.525 mm at each end. A Kepco 0-100 V, 0-5 A power supply provided DC power to 
the heater. A 250 Watt, 2±0.01 Q precision resistor was connected in series with the 
heater. The purpose of the precision resistor was to determine the power to the heater. 
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A schematic of the power distribution system is given in Figure 10. The bottom 
of the base plate was mounted on foam rubber insulation. The sides and the top surface 
of the plexiglas cover were also covered with insulation to reduce the heat losses to 
ambient. 



The voltage drops across the precicision resistor and the heater were measured 
by the Data Acquisition/Control Unit. Since the heater and the precision resistor were 
connected in series the electrical current passing through the resistor was equal to the 
current through the heater. The power to the heater was determined using the following 
relationships, where V/, , vm the measured DC voltage across the heater and V/, , the 
voltage drop across the precision resistor with resistance R. 
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Equation 2.2 


Power = I hxVh Equation 2.3 

2. Thermocouples 

The temperature of the aluminum base plate was measured by using 23 Copper 
Constantan (T-Type) cement-on thermocouples. They were placed on the bottom surface 
of the finned base plate as seen in Figure 11. Since the thermocouples were placed 
between the heater pad and the base plate, flat ribbon thermocouples were chosen to 
provide better surface contact. The thermocouples were glued using Omegabond 101 
which was a high thermal conductivity epoxy. 

Eight Copper Constantan (T-Type) thermocouples were fed through the holes 
drilled into the plexiglas cover to measure the inlet and outlet coolant temperatures. Two 
thermocouples were affixed between the plexiglas cover and the insulation layer w^ere 
the two others were mounted on top of the insulation layer to be able to estimate the heat 
losses through the cover and the insulation. 

Another Copper Constantan (T-Type) thermocouple was placed into the 
refrigerated bath to measure the coolant temperature in the bath. All the thermocouples 
were then connected to the Data Acquisition System. 

3. Data Acquisition System 

All the thermocouples and voltages were measined by an Hewlett Packard 3852 
Data Acquisition/Control Unit. Two 24-channel, thermocouple compensated high speed 
multiplexed boards (HP 44713 A) measured a total number of 34 thermocouple jimction 
temperatures and a 24-charmel high speed multiplexed board (HP 44705 A) measured the 
voltage drops across the flowmeter, precision resistor and the power supply. 





figure 11. Thermocouple Implementation. 
























All theimocouple extensions were labeled and plugged into channel numbers 0 
through 17 of both HP 44713 A boards and the voltage drops were measured by the 
channels 0 through 2 of the HP 44705 A board. 

The Data Acquisition Unit was controlled by a Personal Computer. An interface 
board, G PIB-PCIIA IEEE-488, manufactured by National Instruments was installed into 
this computer. A program attached as Appendix B was written in QBASIC code to store 
the data for later computation and display. A schematic of the entire experimental 
apparatus is presented in Figure 12. 
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Figure 12. The Entire Experimental Apparatus. 


1- Refrigerated Tank 

2- Pump 

3- Pump Speed Control 

4- Flowmeter 

5- Manifold 

6- Heat Exchanger 

7- Power Supply 


8- Precision Resistor 

9- Data Acquisition Unit 

10- Thermocouples 

11- Personel Computer 
VO- Voltage (Flowmeter) 
VI-Voltage (Heater) 

V2- Voltage ^.resistor) 
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DL EXPERIMENTAL PROCEDURE 


A. HEAT TRANSFER EXPERIMENTS 

Heat transfer experiments were performed for both water and Brayco Micronic 
889 (Castrol) Polyalphaolefm (PAO). Since the purpose of these experiments was to 
determine the Colburn j factor for varying Reynolds numbers, the flowrate was varied to 
be able to obtain the desired Reynolds numbers of interest. 

The flow velocity (Reynolds number) was limited by the maximum pressure the 
heat exchanger could withstand. With the power to the heater held constant, the pump 
speed was varied between 10-70 % of its maximum value which corresponded to 
Reynolds numbers in the range of 50-650 for water and Reynolds numbers in the range of 
10-80 for Brayco Micronic 889 Polyalphaolefm. This procedure was repeated for various 
power inputs to the heater (100-250 Watts). The power input to the heater was limited to 
250 Watts due to precision resistor limitations. 

The bath temperature was maintained at around 20‘’C. As stated previously, 
after the system reached the steady stale condition, the necessary data was collected by 
the E>ata Acquisition Unit and used to determine the heat transfer characteristics of the 
test section. For the calculations the following steps were followed: 

1. Reynolds Number Calculations 

All the calculations were based on the hydraulic diameter. The hydraulic 
diameter relates the cross-sectional area of a channel to its heat transfer area. For the 
present study different area deflnitions were defined to determine the hydraulic diameter 
and the Reynolds number. 

• Chaimel Cross Sectional Area (/4 c): 

The product of the fin height h, and the fin spacing s, gives the cross-sectional 
area of each channel. 
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Ac = sy.h 


Equation 3.1 


• Total Cross Sectional Area (A /): 

The total cross sectional area of the test section can be obtained by multiplying 
the number of channels M in the array with Ac. There are nineteen and a half channels 
along the first fin row. 

A/ = MxAc Equation 3.2 

• Heat Transfer Area {A ); 

The heat transfer area of a channel can be determined by using the following 
formula where /i is the fin height, / is the fin length and / is the fin thickness. 

^ = 2 X (A X 0 + 2 X 0 X h)+s X / Equation 3.3 

Since the plexiglas cover is not considered a heat transfer surface, the channel 
length was modeled in Equation 3.3 as an uncovered channel, and included the blunt 
edges of the alternating fins. 

• Test Section Heated Area (Ah): 

The effective area of the heater was defined as the heated area which is actually 
0.07742 m=. 

For the present study the hydraulic diameter was defined as follows; 

4/4 c 

Dh = Equation 3.4 
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All represents the characteristics length of each channel. As stated previously, 
the volumetric flow rate, Q was determined by using the flowmeter calibration curves 
obtained for both water and Brayco Micronic 889 Polyalphaolefin. (See Appendix A). 
Dividing the volumetric flow rate by the total cross-sectional area A /, the velocity V 
through the channels was determined. 


Equation 3.5 


The fluid properties were determined by using the film temperature, 7/. The 
density p, viscosity p, specific heat Cp, thermal conductivity kf and Prandtl number Pr 
of water taken from Ref 9 were curve fitted in the range of 0-80 "C.The same procedure 
was repeated for Brayco Micronic 889 Polyalphaolefin. The properties of this coolant are 
tabulated as Appendix C according to the data provided by the manufacturer, Castrol 
[Ref. 10]. 

Finally the Reynolds number based on the hydraulic diameter was determined 
as follows; 


Rep 




Equation 3.6 


2. Colburn j Factor Calculations 

The film temperature Ti , was averaged between the plate average temperature 
Tavg, and the coolant inlet temperature Ti„. As explained in the section about 
thermocouples, the inlet temperature of the coolant was measured at four locations. The 
plate average temperature was the average temperature of the 23 thermocouples which 
were placed between the heater pad and the base plate. 


Tf 


2 


Equation 3.7 


Tp, 


A15 


Equation 3.8 
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?"(ioo-iii)+7’(40(mio) 
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Equation 3.9 


(The subscript numbers of T in Equations 3.8 and 3.9 represent the thermocouple 
numbers). 

The electrical power input PtUc to the heater was calculated by using the 
equations 2.1, 2.2 and 2.3. The temperature difference through the insulation layer AT ^„,, 
used to estimate the heat losses from the heat exchanger. 


Ar™,,. = 


T416+T417-Ti16-7’i17 

2 


Equation 3.10 


Since the thickness of the insulation tmsui was about 0.0190S m and the thermal 
conductivity kj^i was 0.04 W/mK. The lost or gained heat was determined by using the 
following formula. 


^ kh,sui>‘^^Ti„suiMh 
^msul. 


Equation 3.11 


Here A/, represents the area covered with the insulation layer. The actual power 
P, to the system was obtained by adding AP to the electrical power. 


P = PeUc +AP Equation 3.12 

The next step was to calculate the heat flux. Since the heater dimensions and 
the power input to the heater was held constant for a set of experiments, the heat flux q” 
was also constant. The heat flux and the average Nusselt number Nuavg were calculated 
as follows. 

p 

q”=~ Equation 3.13 


26 



_ 


Equation 3.14 


Finally, the experunental Colburn j factor was calculated to relate the 
dimensionless heat transfer coefficient to the Reynolds number 


’ ReoxPr'^ 


Equation 3.15 


B. FLOW VISUALIZATION EXPERIMENTS 

The procedure for these experiments was to change the flow velocity of the 
water and visualize the flow on the offset strip fins. Only water was used as the working 
fluid and the power input to the heater was set to zero during the visualization 
experiments. The pump speed was varied to obtain the same Reynolds numbers of the 
heat transfer experiments performed with water (Re 50-650). Each fin row was labeled 
on the plexiglass cover to be able to determine the location of the changes in the velocity 
profile. First blue ink and then white dye were injected into the plexiglas cover and 
directed through the fin array. Photographs were taken for each Reynolds number 
settings. Several photographs should be taken to be able to observe all the fin rows. 
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IV. RESULTS 


A. HEAT TRANSFER RESULTS 

When water was used as the working coolant the power supply was set to 100, 
150, 200 and 250 Watts. The data obtained and the calculated Colburn j factors are 
tabulated for each power setting in Tables 1-4. Also some of the experimental data were 
plotted for varying Reynolds numbers and power inputs. 

The differences between the water inlet, and water outlet temperatures are 
plotted on a graph (Figure 13), while the differences between the average plate and water 
inlet temperatures are plotted on another graph (Figure 14), for all power settings. Figure 
15, shows the plate average temperatures for varying Reynolds numbers. The calculated 
Colburn j factors are also plotted separately for each power setting in Figures 16-19. 

The same procedure was repeated for the coolant Brayco Micronic 889 
Polyalphaolefin. This time the power supply was set to 100,175 and 250 Watts. The data 
obtained were tabulated in Tables 5-7. Temperature differences mentioned above are 
plotted in Figures 20 and 21. The calculated Colburn j factors for the coolant PAO are 
shown in Figures 22-24. The correlations developed by Weiting [Ref 3], and Manglik 
and Bergles [Ref 4], were used to predict the Colburn j factor for the present study. (See 
Chapter 1, Section C). These predicted Colburn j factors were plotted in Figures 25-28 to 
be able to compare with the experimental results obtained for various power inputs. 
These figures also included the results for the coolant PAO to determine the Prandtl 
number effect on the Colburn j factor. It was noted that the Prandtl number of the PAO 
(Brayco Micronic 889) was aroimd 230-250 and the Prandtl number of water was around 
5-7 during the experiments. 

The data and the graphs showed that as the Reynolds number increased, the 
Nusselt number also increased for both of the coolants. It was noted that higher power 
inputs produced higher Nusselt numbers. For the same power input and Reynolds 
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number, the coolant Brayco Micronic 889, produced significantly higher Nusselt 
numbers than water. 

As the Reynolds number increased, the temperature difference between the inlet 
and outlet of the heat exchanger showed an exponential decrease for all power settings. 
Since the inlet temperature of the coolant was maintained at a constant temperature 
(around 20 ®C), for the same Reynolds number as the power input increased the 
temperature difference mentioned above also increased. 

As the power input increased, the average surface temperature of the base plate 
increased. But increasing the Reynolds number decreased the plate average temperature 
and also decreased the difference between the plate average and coolant inlet 
temperatures. Since the volume of the circulating bath was small, as the Reynolds 
number exceeded about 650 for water and about 60 for the PAO, the inlet temperature of 
the coolant and the plate average temperature began to increase. When the coolant was 
the PAO, for the Reynolds numbers greater than 60, the maximum increase in the plate 
average temperature was less than 0.60 °C. 

It was noted that, even though the variation was very small, higher power inputs 
produced higher Colburn j factors for the same Reynolds number settings. As the power 
input increased the plate average temperature and the coolant film temperature also 
increased. This caused a decrease in the Prandtl number and an increase in the Colburn j 
factor. The predicted Colburn j factors developed by Weiting [Ref 3], and Manglik and 
Bergles [Ref 4], were compared to the experimental results. When the water was used as 
the working fluid, there was a very close agreement with the results of these studies. 
Especially, when the Reynolds munber was less than 200, the experimental results were 
in the accuracy zone of Weiting (±10 %), and Manglik and Bergles (± 20 %). For the 
Reynolds numbers greater than 200, the experimental results deviated slightly from the 
accuracy zone. The maximum deviation was 38.20 %, when the power input was 
minimum (100 Watts), and 29.93 %, when the power input was maximum (250 Watts). 
For the power input of 250 Watts, about half of the experimental results were in the 
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accuracy zone. The difference between the predicted and calculated Colburn j factors for 
the power input of250 Watts were tabulated in Table 8. 

Since the predictive correlations were applicable for Prandtl numbers in the 
range of 0.5 to 15, and the Prandtl number of Brayco Micronic 889 (PAO), was about 
230-250 during the experiments, there was not a close agreement between the predicted 
and experimental results when this coolant (PAO) was used as the working fluid. Even 
though the experimental results were higher than the predicted values the curves obtained 
for both of the experimental and predicted Colburn j factors showed a similar tendency. 
(See Figures 25-28). This validated the Prandtl number limitations of these predicted 
correlations. 

Finally, the Colburn j factors of water and Brayco Micronic 889 (PAO) were 
compared to each other. When the power input was 250 Watts and the Reynolds number 
was around 57, the calculated Colburn j factor was 0.0856 for water (Pr = 6.03), and 
0.3280 for the PAO (Pr = 238.72), it was noted that the ratio between the Colburn j 
factors was about the same to the 1/3 to the power of the ratio between the Prandtl 
numbers. 




The difference between the ratios was about 11 %. 
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Powey Input=lBa Watts <Heat Flux g1299 > 


Re 

T1 

(hhtTtnp.) 

Cci 

T2 

(*C) 

T3 

(AofPM*) 

rc) 




g 

g 

S3. 01 

20.22 

20.80 

25.03 

QB 

IBI 




142.19 

20.35 

20.74 

24.91 

OBI 

bob 

BB 

BIB 

IBI 

186.11 

20.28 

20.62 

24.70 

EBI 

^^B 




229.61 

20.25 

20.55 

24.42 

OBI 

^S9 



IBI 

273.12 

20.18 

20.43 

24.28 

BBI 

BBI 

BBI 

bb 

bbb 

316.65 

20.05 

20.26 

24.02 


I^BI 


bei 

bbi 

359.23 

20.06 

20.24 

23.94 

BB 

^^B 

HB 

bb 

QIBI 

402.83 

20.09 

20.23 

23.86 


jj^BI 

IBI 


B^l 

447.02 

20.19 

20.30 

23.82 

BQI 

BBI 

^^^1 

BH 

B^l 

490.96 

20.26 

20.36 

23.76 

mgm 


0.560 

10.13 

^bb 





0.09 

BBI 




BBS 


BB 

BBI 


Bl 

I^BB 

I^B 



Table 1. Experimental Data for Water (Power Input = 100 W). 








Powe^Injgut^jl50Jlatt8<Heat^lin<®^94^Uj2|^ 


Re 

T1 

(WMTtMp.) 

Co 

T2 

[OmMTmip) 

CC) 

t3 

rc) 

T2-T1 

Ct) 

T3-T1 

rc) 

(Watts) 

HV 

J 

S 4.21 

20.11 

20.96 

26.89 

0.85 

6.78 

0.527 

7.83 

0.0779 

150.40 

20.00 

20.59 

26.11 

0.59 

6.11 

0.592 

8.69 

0.0325 

188.58 

20.02 

20.52 

25.98 

0.50 

5.96 

0.605 

8.93 

0.0255 

: 232.81 

19.95 

20.39 

25.79 

0.44 

5.84 

0.611 

9.11 

0.0210 

ZTi.ZZ 

19.95 

20.33 

25.69 

0.38 

5.74 

0.619 

9. 28 

0.0180 

321.81 

19.98 

20.31 

25.58 

0.33 

5.60 

0.629 

9.51 

0.0159 

36 S .99 

19.94 

20.22 

25.50 

0.28 

5.56 

0.074 

9.60 

0.0141 : 

410.64 

19.97 

20.24 

25.46 

0.26 

5.49 

0.631 

9.72 

0.0127 

; 455.44 

20.06 

20.29 

25.39 

0.23 

5.33 

0.631 

10.00 

0.0118 

498.72 

20.06 

20.26 

25.17 

0.20 

5.11 

0.629 

10.45 

0.0112 

543.72 

20.17 

20.35 

25.12 

0.18 

4.95 

0.631 

10.77 

0.0106 

58*'. 64 

20.17 

20.34 

25.03 

0.17 

4.86 

0.632 

10.96 

0.0100 


Table 2. Experimental Data for Water (Power Input =150 W). 

























Powei* Input =200 Watts <Heat Flux - 2580 


Re 

T1 

t*C) 

T2 

Cc) 

T3 

(Av9Ptot*) 

Cc\ 

T2-T1 

("C) 

T3-T1 

rn 

(Watts) 

Hit 

3 

56.03 

20.50 

21.66 

29.10 

1.16 

8.60 

0. 390 

8.16 

0.0795 

149.19 

20.67 

21.50 

28.35 

0.83 

7.68 

0.475 

9.15 

0.0334 

195.36 

20.63 

21.36 

28.14 

0.73 

7.51 

0.478 

9.38 

0.0261 

241.66 

20.63 

21.27 

28.03 

0.64 

7.40 

0.483 

9.52 

0.0214 

28''.S4 

20.58 

21.14 

27.89 

0. 56 

7.31 

0.484 

9.63 

0.0182 

333.41 

20.61 

21.09 

27.72 

0.48 

7.11 

0.493 

9.92 

0.0161 

379.60 

20.68 

21.12 

27.61 

0.44 

6.93 

0.494 

10.17 

0.0146 

425.53 

20.71 

21.10 

27.50 

0.39 

6.79 

0.494 

10.39 

0.0132 

472.07 

20.81 

21.17 

27.44 

0.36 

6.63 

0.503 

10.64 

0.0122 

518.14 

20.88 

21.24 

27.33 

0.33 

6.45 

0.507 

10.91 

0.0114 

562.16 

20.82 

21.14 

27.07 

0.32 

6.25 

0.515 

11.27 

0.0109 

606.96 

20.86 

21.17 

26.86 

0.31 

6.00 

0.519 

11.74 

0.0105 


Table 3. E^qjerimental Data for Water (Power Input = 200 W). 
























Power 


nput=2SB Watts <Heat Flux « 323B W/w^> 


Re 

T1 

(hktTMip.) 

T2 

(CMhtTtap) 

Cn 

T3 

t*C1 

T2-T1 

("C) 

T3-T1 

rc) 


Mu 

3 

S'?. 04 

20.59 

22.01 

30.44 

1.42 

9.85 

0.543 

8.89 

0.0856 

152.41 

20.65 

21.68 

30.07 

1.03 

9.42 

0.548 

9.31 

0.0335 

199.61 

20.62 

21.53 

29.86 

0.91 

9.24 

0.554 

9.50 

0.0261 

246.24 

20.56 

21.39 

29.59 

0.83 

9.03 

0.564 

9.71 

0.0216 

292.56 

20.57 

21.31 

29.27 

0.74 

8.70 

0.569 

10.09 

0.0189 

338.04 

20.49 

21.16 

28.93 

0.67 

8.44 

0.577 

10.41 

0.0168 

384.29 

20.50 

21.09 

28.76 

0.59 

8.26 

0.590 

10.65 

0.0151 

430.94 

20.60 

21.13 

28.61 

0.53 

8.01 

0.602 

10.99 

0.0139 

47'?. 48 

20.64 

21.11 

28.51 

0.47 

7.87 

0.615 

11.16 

0.0127 

523.55 

20.72 

21.14 

28.31 

0.42 

7.59 

0.632 

11.58 

0.0121 

570.20 

20.81 

21.20 

28.20 

0.39 

7.39 

0.642 

11.91 

0.0114 

614.56 

20.74 

21.12 

27.96 

0.38 

7.22 

0.645 

12.20 

0.0108 


Table 4. Ejqperimental Data for Water (Power Liput = 250 W). 
























Power Input-lBfl Watts <Heat Flux ^1290 W/n^> 


£e 

T1 

((ftMTMp.] 

C°cj 

T2 

rcj 

T3 

(A««I>M«) 

<*C) 

T2-T1 

rc) 

T3-T1 

rc) 


Hu 

3 

13.48 

20.49 

23.36 

26.17 

2.87 

5.68 


85.49 




21.96 

25.35 

1.89 

5.28 






21.28 

24.87 

1.36 

4.95 








1.01 

4.62 



1^33 





0.75 

4.40 

0.665 

110.64 

0^3 




3IQ 

0.62 

4.26 


QQI 




0^31 


0.60 

4.16 

0^31 

Q^3 






0.53 

4.01 


QQ 

^^3 

L^n 

IKI 

Bffii 

ESI 

0.52 

3.98 


E^B 



Table 5. Experimental Data for Brayco Micronic 889 (Power Input =100 W). 












Power Input-175 Watts <Heat Flux « 2260 W/w»> 


He 

T1 

(*C) 

Ti 

T3 

(Ay9J>tM4) 

T2-Tj. 

iX) 

T3-T1 

(’c» 

(WHts) 

■u 


13.69 

20.28 

2S.13 

29.56 

4.85 

9.28 

-0.124 

91.33 

1.0675 

22. IS 

19.92 

23.32 

28.68 

3.40 

8.76 

0.078 

96.85 

0.6982 

31.08 

19.79 

22.13 

27.74 

2.34 

7.95 

0.171 

106.76 

0.5476 

39.88 

19.67 

21.55 

27.20 

1.88 

7.53 

0.242 

112.75 

0.4502 

48.12 

19.7S 

21.18 

26.92 

1.43 

7.17 

0.327 

118.46 

0.3920 

SS.38 

19.86 

21.05 

26.87 

1.19 

7.01 

0.347 

121.18 

0.3484 

61.37 

19.93 

20.98 

26.84 

1.05 

6.91 

0.376 

122.96 

0.3190 

69. 26 

20.32 

21.27 

27.07 

0.95 

6.75 

0.386 

125.89 

0.2897 

72.13 

20.73 

21.64 

27.37 

0.91 

6.64 

0.378 

127.98 

0.2831 


Table 6. Eaiperimental Data for Brayco Micronic 889 (Power Input = 175 W). 





















Power Input=250 Watts <Heat Flux ^ 3230 W/ra*> 


Re 

T1 

(•c) 

T2 

C"C) 

T3 

(A«9Pltl«) 

(*C) 

T2-T1 

rcj 

T3-T1 

(•cj 



B 

13.70 

20. 54 

27.64 

33.76 

7.10 

13.22 




22.73 

20.96 

25.61 

33.01 

4.65 

12.05 








IQQ 


Q^l 
















2.16 









1.86 









1.72 












1^20 


B9 

EQ 

JESS 

B3 

IQ 

I^B 


BD 



Table 7. Experimental Data for Brayco Moronic 889 (Power Input = 250 W). 
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Figure 13. Temperature Difference vs. Reynolds Number. 







Figure 14. Temperature Difference vs. Re}molds Number 
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Figure 15. Plate Average Temperature vs. Reynolds Number. 















Figure 17. Colburn j Factor vs. Reynolds Number (Power = 150 "W). 









figure 18. Colburn j Factor vs, Reynolds Number (Power = 200^). 




Figure 19. Colburn j Factor vs. Reynolds Number (Power = 250 "W). 










Figure 20. Temperature Difference vs. Re 3 molds Number. 
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Figure 22. Colburn j Factor vs. Reynolds Number (Power = lOD W). 
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Figure 23. Colburn j Factor vs. Reynolds Number (Power = 175 W). 










figure 24. Colburn j Factor vs. Reynolds Number (Power = 250"W). 
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Figure 25. Comparison of Predicted and Experimental j's. 
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Comparison of Predicted and Experimental j's. 






f le^uamusdx^ pue uosfreduzoQ -g^ 
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Powttr Input-25B Watts 







RbymoXds 

MuMbar 

3 

i Dx€£erance 
(Wdibxng) 

( Dx£facanca 
(Mangiik 4 Becglaa] 

152.41 

0.0335 

-8.8705 

-2.4176 

1S9.61 

0.0261 

2.0190 

7.7195 

246.24 

0.0216 

9.3036 

14.5032 

252.56 

0.0189 

13.1094 

18.0300 

: 338.04 

0.0168 

16.2434 

20.9368 

: 384.29 

0.0151 

19.3857 

23.8617 

430.94 

0.0139 

21.1444 

25.4860 

477.48 

0.0127 

23.7010 

27.8701 

523.55 

0.0121 

24.1215 

28.2391 

570.20 

0.0114 

25.0518 

29.0929 

614.56 

0.0108 

25.9607 

29. 9303 


Table 8. ^ Difference Be-tween The Predicted and 
Calculated j Factors. (Power = 250 W). 

























B. FLOW VISUALIZATION RESULTS 

In this study, a large number of photographs were taken to be able to visualize 
the flow through the offset strip fins. Two of them were presented in Figures 29 and 30. 

When the ink streak was introduced at the leading edge of the fin, the ink 
followed the side of the fin and did not rise perceptibly firom the level of injection. The 
ink then curled to its center, around the trailing edge of the fin and followed the open 
channel until reaching the leading edge of the next row. Afier being deflected by the 
laminar boundary layer forming on the fin, the flow pattern repeated itself (See Figures 
29 and 30). Since the Reynolds number was limited to 650 by the maximum pressure the 
heat exchanger could withstand, transition to turbulent and vortices were not observed. 
But, it was noted that as the Reynolds number increased, the streaklines produced a more 
severe wavy pattern. 



Figure 29. Flow Visualization of Water (Re = 200, First Five Rows). 
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figure 30. Flow Visualization of Water (Re = 200, After First Five Rows). 
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V. CONCLUSIONS 


The heat transfer experiments of this study showed that, when the working fluid 
was water, the experimental values come to close agreement to the predicted values. But 
when Brayco Micronic 889 (PAO) was used as the working fluid, the experimental 
results were significantly higher than the predicted results. This provides evidence that 
the previous correlations were valid for water but not for PAO. 

The comparison of the predicted dimensionless heat transfer coefficients to the 
experimental values validated the Prandtl number dependency of the Colburn j factor. It 
was obtained that the increase in the Colburn j factor was about the same with 1/3 to the 
power of the increase in the Prandtl number of the working fluid. This also proved the 
superiority of the higher Prandtl number fluids on the lower Prandtl number fluids to 
increase the performance of an heat exchanger. 

Flow visualization experiments showed that, the interrupted geometry of the 
fins prevented the velocity boundary layers firom getting excessively thick. Observing 
these laminar boundary layers along the offset strip fins helped to explain the 
enhancement in the heat transfer of the system for varying flow rates. 

In general, this study will help the researchers to predict the dimensionless heat 
transfer coefficient of a ten times scaled up model of a liquid flow through module used 
in electronics cooling. 
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VL RECOMMENDATIONS 


In any continuation of this study the followings are reconunended; 

• Clean the water deposits on the base plate and offset strip fins and prepare the 
test section for the new study. 

• Redesign and manufacture a stronger, pressure resistant plexiglas cover for the 
heat exchanger. 

• Insert additional thermocouples into the inlet and outlet of the test section. 
Remove the heater and calibrate the thermocouples mounted on the bottom of the base 
plate fin assembly to a zero reference point. 

• Repeat the experiments of the present study by using a pump with higher 
capacity to investigate the heat transfer and flow characteristics of the system for the 
turbulent flow. 

• Develop a numerical solution to the flow through the offset strip fins and 
compare the solution with the experimental results. 

• Investigate the Fanning friction factor and the effect of the fin thickness on the 
form drag. 
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APPENDIX A. FLOWMETER CALIBRATION 


The flowmeter was calibrated in two different ways. The first one consisted of a 
comparison of the average voltage output to the average time to collect a specific 
volume of water as shown in Figure 31. 



Fiflure 31. Flowmeter Calibration (Coolant: Water). 
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The second method was based on the collected mass of the water in a specific 
time for varying pump speeds as shown in Figure 32. 


Fknmeter Curve FItiirg 



0 20 « 60 80 10D 


I _ FIow(%) _ I 

Figure 32. Flovnneter Calibration (Coolant: Water). 
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The calibration of the flowmeter for Brayco Micronic 889 Polyalphaolefin (PAO) 
was based on the first method which consisted the comparison of the flow percentage 
readings to the average time to collect a specific volume of PAO as shown in Figure 33. 



Piaure 33. Flowmeter Calibration (Coolant; PAO), 
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APPENDIX B. QBASIC PROGRAM CODE 


CLS 

COLOR 14,9 

PRINT" Program for National Instruments PC2 IEEE-488 BOARD TO HP-3852" 
PRINT" This uses the Universal Language Interface" 

PRINT " ULI.COM must be run prior to running the program" 

PRINT" The program name is COLBURNJ.BAS” 

PRINT" Modified by Ata KOMRAL" 

OPEN "GPIBO" FOR OUTPUT AS #1 
OPEN "GPIBO" FOR INPUT AS #2 
PRINT #1, "ABORT" 

PRINT #1, "RESET" 

PRINT # 1, "GPIBEOS CR LF" 

PRINT #1 "CLEAR 9" STARTIT' 

DIM TEMP(20). TMP(20), V(3), TIN(33), TOUT(33), DIFF(33), TAV(33), FILM(33), 
TINSULA(33). TPGLSS(33) 

PRINT #1, "REMOTE 9" 

PRINT #1. "OUTPUT 9;RESET 700" 

PRINT #1, "OUTPUT 9;USE 700" 

INPUT "ENTER PUMP SETTING ", SET 
TOTAV = 0! 

TOTIN = 0! 

TOTOUT = 0! 

TOTDIFF = 0! 

TOTFTLM = 0! 

TOTCVR = 0! 

TOTINSL = 0! 

N = 33 

PRINT TIMES; 

FORJ=lTON 

PRINT #1, "OUTPUT 9;CONF TEMPT " 

FORI = 0TO17 
CHNL$ = STR$(I+ 100) 

MSGS = "OUTPUT 9;MEAS TEMPT" + CHNLS 
PRINT #1, MSGS 
PRINT #1, "ENTER 9" 

INPUT #2, TMPS 
TMPa) = VAL{TMP$) 

NEXTI 

PRINT #1, "OUTPUT 9; CONF TEMPT" 

FORI = 0TO17 
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CHN$ = STR$(I + 400) 

MSGES = "OUTPUT 9; MEAS TEMPT" + CHN$ 
PRINT #1,MSGE$ 

PRINT #1, "ENTER 9" 

INPUT #2, TEMPS 
TEMP(I) = VAL(TEMP$) 

NEXTI 

PRINT #1, "OUTPUT 9;CONF DCV" 

PRINT #1, "OUTPUT 9;USE 700" 

FOR I = 0 TO 2 
CHS = STRS{I + 500) 

MESSS = "OUTPUT 9;MEAS DCV" + CHS 
PRINT #1,MESSS 
PRINT #1, "ENTER 9" 

INPUT #2, VOLTS 
V(I) = VAL(VOLTS) 

NEXTI 

PRINT #1, "CLEAR 9" 

PRINT #1, "LOCAL 9" 

TPLATEl = 0! 

1 = 0 

DO WHILEI< 12 
1 = 1 + 1 

TPLATEl = TPLATEl + TMP(I - 1) 

LOOP 

TPLATE2 = 0! 

1 = 0 

DOWHILEI< 11 
1 = 1+1 

TPLATE2 = TPLATE2 + TEMP(I - 1) 

LOOP 

TSURF = (TPLATEl + TPLATE2) / (12 + 11) 
TIN = 0! 

1 = 0 

DOWHILEI<4 
1 = 1 + 1 

TIN = TIN + TEMP(I+ 11) 

LOOP 

INTMP = TIN/4 
TOUT = 0! 

1=0 

DOWHILEI<4 
1 = 1 + 1 
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TOUT = TOUT + TMP(I + 11) 

LOOP 

OUTTMP = TOUT/4 
TPGLASS = 0! 

1 = 0 

DOWHILEI<2 
1 = 1 + 1 

TPGLASS = TPGLASS + TMP(1 + 15) 
LOOP 

TCVR = TPGLASS/2 
TINS = 0! 

1 = 0 

DOWHILEI<2 
1 = 1+1 

TINS = TINS + TEMPO + 15) 

LOOP 

TINSUL = TINS/2 
TAV(J) = TSURF 
TIN(J) = INTMP 
TOUT(J) = OUTrMP 
DIFF(J) = OUTTMP - INTMP 
FILM(J) = (TSURF + INTMP) / 2 
TPGLSS<J) = TCVR 
TINSULA(J) = TINSUL 
TOTAV = TOTAV + TAV(J) 

TOTIN = TOTIN + TIN(J) 

TOTOUT = TOTOUT + TOUT(J) 
TOTDFF = TOTDIFF + DIFF(J) 

TOTFILM = TOTFILM + FILM(J) 
TOTCVR = TOTCVR + TPGLSS(J) 
TOTINSL = TOTINSL + TINSULA(J) 
NEXTJ 

TSURFACE = TOTAV / N 

INLET = TOTIN/N 

OUTLET = TOTOUT/N 

DELTA = TOTDIFF/N 

FILMTEMP = TOTFILM /N 

TCOVER = TOTCVR / N 

TINSULATION = TOTINSL / N 

DELTINSUL = TINSULATION - TCOVER 

TEMDIF = TSURFACE - INLET 

R = 2 

P = V(2)*(V(1)/R) 
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FOR I = 0 TO 2 

PRINT USING "###.##"; V(I); 

NEXT I 

PRINT "PWR=" 

PRINT USING "##.##"; P 
FORI = 0TO3 
PRINT I + 100; 

PRINT USING "###.###"; TMT®; 

PRINT I + 104; 

PRINT USING "###.###"; TMP(I + 4); 

PRINT I + 108; 

PRINT USING "###.###"; TMP(I + 8); 

PRINT I + 112; 

PRINT USING "###.###"; TMP(I + 12) 

NEXTI 
FORI = 0TO3 
PRINT I + 400; 

PRINT USING "###.###"; TEMPa); 

PRINT I + 404; 

PRINT USING "###.###"; TEMP(I + 4); 

PRINT I+ 408; 

PRINT USING "###.###"; TEMPa + 8); 

PRINT I+ 412; 

PRINT USING "###.###"; TEMP(1 + 12) 

NEXTI 

PRINT # 1, "OUTPUT 93EEP" 

PRINT #1, "CLEAR 9" 

PRINT #1, "LOCAL 9" 

z= 19.5 

1 = .03175 

s=.01175 

t=.00152 

h=.01524 

Ah = .07742 

Ac = s * h 

Af = Ac * z 

At = s*l + 2*h*l + 2*t*h 
Dh = (4 * Ac) / (At /1) 

vis = -{3.7292 /10 12) * FILMTEMP 3 + (7.125 / 10 10) • FILMTEMP 2 
- (5.0808 / 10 ^ 8) * FILMTEMP + (1.7758 /10 ^ 6) 

kf = (8.3333 / 10 8) * FILMTEMP ^ 3 - (2.1429 /10'' 5) * FILMTEMP ^ 2 
+ (2.631 /10 3) • FILMTEMP + .55206 

Pr = (4.4271 / 10 ^ 7) * FILMTEMP ^ 4 - (1.0604 /104) * FILMTEMP ^ 3 
+ (9.9979 / 10 ^ 3) * FILMTEMP ^ 2 - (4.9008 /10) ♦ FILMTEMP + 13.6 
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den = (2 /10 5) * FILMTEMP 3 - (6.407 /10 ^ 3) * FILMTEMP 2 
+ (3.207 /10 ^ 2) * FILMTEMP + 1002.3 

Cp = (5.2604 /10'' 6) * FILMTEMP ^ 4 - (1.1167 /10 ^ 3) * FILMTEMP 3 

+ (9.3021 /10 2) • FILMTEMP'' 2 - (3.2558) * FILMTEMP + 4217.8 

m = .0018* SET- .0073 

Q = ni/den 

VEL = Q/Af 

Re = (VEL*Dh)/vis 

kins = .04 

thick = .75 • 2.54 /100 

DELP = kins * Ah • DELTINSUL / thick 

PTOT = P + DELP 

FLUX = P/Ah 

Nu = (FLUX • Dh) / kf / TEMDIF 
Je = Nu/Re/(Pr^(l/3)) 

PRINT TIMES 
PRINT '’TAVG="; 

PRINT USING "###.###"; TSURFACE; 

PRINT "TIN="; 

PRINT USING "###.###"; INLET; 

PRINT "TOUT="; 

PRINT USING "###.###"; OUTLET; 

PRINT "DELTA="; 

PRINT USING "###.###"; DELTA; 

PRINT "TFrLM="; 

PRINT USING" ###.###"; FILMTEMP; 

PRINT "TDIF="; 

PRINT USING" ###.###"; TEMDIF; 

PRINT "TCOVER="; 

PRINT USING " ###.###"; TCOVER; 

PRINT "TINSUL="; 

PRINT USING" ###.###"; TINSULATION; 

PRINT "visc=''; 

PRINT USING "###.######"; vis; 

PRINT "cond="; 

PRINT USING "##.######"; kf; 

PRINT "Prandtl="; 

PRINT USING "##.######"; Pr, 

PRINT "Density="; 

PRINT USING "####.#####"; den; 

PRINT "Cp="; 

PRINT USING "####.####"; Cp; 

PRINT ’T>OWER="; 
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PRINT USING "####.####'•; PTOT; 

PRINT "m="; 

PRINT USING "###.####"; m; 

PRINT "VEL="; 

PRINT USING "###.#####"; VEL; 

PRINT "Re="; 

PRINT USING "###.#####"; Re; 

PRINT "FLUX="; 

PRINT USING "####.#####"; FLUX; 

PRINT "Nu="; 

PRINT USING "####.#####"; Nu; 

PRINT "j="; 

PRINT USING "####.#####"; Je; 

INPUT " DO YOU WISH TO SAVE DATA ? ", SAVITS 
IF SAVITS = "n" THEN GOTO QUES 
IF SAVITS = "N" THEN GOTO QUES 
FILEIT: 

PRINT " IF DATA FILE EXISTS OLD DATA WILL BE OVERWRITTEN " 
INPUT "ENTER THE NAME OF THE DATA FILE TO CREATE", FILES 
IF FILES ="" THEN GOTO FILEIT 
OPENTT; 

OPEN FILES FOR OUTPUT AS #3 
PRINT #3, DATES 
PRINT #3, TIMES 
PRINT #3, "Pump Set ="; 

PRINT #3, SET 
PRINT#3,"EIec.Power="; 

PRINT #3, USING "###.###"; P 
PRINT #3, "Act.Power ="; 

PRINT #3, USING "###.###"; PTOT 
PRINT #3, "Lost ="; 

PRINT #3, USING "###.###"; DELP 
PRINT #3, "Flowrate ="; 

PRINT #3, USING "###.####"; m 
PRINT #3, "Density = "; 

PRINT #3, USING "###.####"; den 
PRINT #3, "Viscosity = "; 

PRINT #3, USING "###.####"; vis 
PRINT #3, "Conductivity = "; 

PRINT #3, USING "###.####"; kf 
PRINT #3, "Prandtl = "; 

PRINT #3, USING "###.####"; Pr 
PRINT#3,"Cp = "; 
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PRINT #3, USING "###.####"; Cp 
PRINT #3, "Velocity ="; 

PRINT #3, USING "###.####"; VEL 
PRINT #3, "Reynolds ="; 

PRINT #3, USING "###.####"; Re 
PRINT #3, "Flux="; 

PRINT #3, USING "####.###"; FLUX 
PRINT #3." Nusselt ="; 

PRINT #3. USING "###.###"; Nu 
PRINT #3, "Colburn j ="; 

PRINT #3, USING "###.####"; Je 
PRINT#3,"Vflmtr="; 

PRINT#3,V(0) 

PRINT#3,"Vheat="; 

PRINT#3.V(1) 

PRINT #3, "Vresis ="; 

PRINT#3.V(2) 

PRINT #3, "Thennocouples(C)" 

FORI = OTO 17 

PRINT #3, (I + 100); TMP(I); (I + 400); TEMPO) 

NEXTI 

PRINT #3, "Tsurf="; 

PRINT #3, USING "###.##"; TSURF 
PRINT#3."Tfilin="; 

PRINT #3, USING "##.###"; FILMTEMP 
PRINT #3, "Tin ="; 

PRINT #3, USING "##.###"; INLET 
PRINT #3. "Tout ="; 

PRINT #3. USING "##.###"; OUTLET 
PRINT #3, "Delta = "; 

PRINT #3, USING "##.###"; DELTA 
PRINT#3,"Tdiff="; 

PRINT #3, USING "##.###"; TEMDIF 
PRINT #3, "Tcover ="; 

PRINT #3. USING "##.###"; TCOVER 
PRINT #3, "Tinsulation ="; 

PRINT #3, USING "##.###"; TINSULATION 
QUES: 

INPUT " DO YOU WISH TO TRY ANOTHER RUN ( Y/N ) ? ", AGAINS 

IF AGAINS = "Y" THEN GOTO STARTIT 

IF AGAINS = "y" THEN GOTO STARTIT 

IF AGAINS = "N" THEN END 

IF AGAINS = "n" THEN END 

GOTO QUES 

END 
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APPENDIX C. PROPERTIES OF BRAYCO MICRONIC 889 (PAO) 


The properties of Brayco Micronic 889 Polyalphaolefin taken from Ref 10 are tabulated 
below. The thermal conductivity of Brayco Micronic 889 (PAO) was 0.083-0.085 BWIhr.ff- 
during the heat transfer experiments. 


TEST(ASTN) 

DESCRIPTION 

RESULT 

D-445 

Kinematic Viscosity 

(10^6 m-'R/s) 


@ 100-C 

1.7 


0 40-c 

52 


0 -40“C 

260 


0 -54 “C 

1200 

D-287 

Density 

(g/ml) 


0 0“C 

0.811 


0 20 *c 

U. 794 


0 40 °C 

0.777 


0 100‘C 

0.723 


0 160“C 

0. 661 

D-2766 

Specific Heat 

(cal/g C) 


0 -18“C 

0.49 


0 10 ®C 

0.52 


0 

0.54 


0 93°C 

0.58 


Table 9. Properties of Brayco Micronic B89[Ref. 10]. 
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APPENDIX D. SAMPLE CALCULATIONS 


The following calculation is for the following setpoint. 

Coolant: Water. 

Electrical Power Input, Peic= 100.172 Watts. 

Flow(%) = 40. 

Mass Flow Rate, m = 0.0647 kg/s. 

Inlet Temperature, T1 = 20.05 ”C 
Outlet Temperature, T2 = 20.26 ®C 
Plate Average Temperature, T3 = 24.02 “C 

Temperature Difference Through the Insulation Layer, AT/mi/ = 3.088 °C 

1. Characteristic Dimension 

Ac = sxh 

Ac = 11.75ffwn X 15.24inm. 

Ac=m.\mm^. 


Af=MxAc 
^/= 19.5x179.1 
Af= 3492.5OTm2. 


A = sxl + 2xhxl+2xhxt 

A = (11.75 X 31.75) + (2 x 15.24 x 31.75) + (2 x 15.24 x 1.52) 
A = 1387.1/«ot^. 


^ 4^ 4x179.1 _ 

^'’"(^//)"(1387.1/31.75) 


2 . 


Properties of Water 


Using Table A.6 [Ref 9]; 
p= 1000.11 kg/m’ 


kf= 0.6005 W/m.K 
Pr = 6.625 


Cp = 4180.51 J/kg.K 
p = 9.623x10-^ N.s/m’ 
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3. Reynolds Number 

m 0.0647 ^ , 

Volumetric Flow Rate, Q = 'p = loOO 11 ^ ^ 


Q 6 469x10“* 

Velocity, V = ^ = = 0.0185 m/s. 

Af 3492.5x10“^ 

pxVxDh (1000.11)x(0.018S)x(16.39x10-3) 
~ 9 623x10-4 


4. Nusselt Number 


Lost/Gained Heat, A/* 


^insul^Ah^^Tj^l 


where kmui (0.04 W/m.K), is the thermal conductivity and /*„„/ (0.01905), is the thickness of the 
insulation layer. 


A/’=- 


-(0.04)x(0.07742)x(3.088) 
0.01905 


•=-0.502 Watts. 


Total Power Input, P = Peic - AP 

P= 100.172 + 0.502= 100.674 Watts. 
P 100.674 


Heat Flux, q'' =— 


Ah 0.07742' 


^ 1300.36 W/m^ 


g”xD^ (1300.36)x(16.39xl0-*) 

” kj^iTi-T\) ~ (0.6005)x(3.97) “ 


Colburn j Factor 

Nuavg 8.94 


/?e/,xPr'^* (315.65)xr6.625)*'^ 
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APPENDIX E. UNCERTAINTY ANALYSIS 


The accuracy of the results of this study was estimated by performing an uncertainty 
analysis. The theorem established by S .J. Kline and F. A. McClintock [Ref 11] was used to 
calculate the uncertainties. In this theorem a calculated result y is a function of several 

independent measured variables, {xi, X 2 ,.,Xn}. Each measured value has some uncertainty, 

{ui, U 2 .,u„}, and these uncertainties lead to an uncertainty in y, which is called u^. To 

estimate Uy, it’s assumed that each uncertainty is small enough that a first-order Taylor 

expansion of y(xi, X 2 ,.x„) provides a reasonable approximation. Under this approximation, 

y is a linear function of the independent variables and the uncertainty is determined by; 



Here, all uncertainties must have the same odds and must be independent of each other. 



The uncertainty analysis performed for the power setting of 100 Watts and the Reynolds 
number of about 315.65 (for water) which produced a Colburn j factor of 0.0151 is presented 
below: 

1. Reynolds Number Uncertainty 
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For T^fo, = 22.04°C 


p = 1000.11 kg/m^ k= 0.6005 W/m.K, n= 9.62 E-4N-s/m^ Pr= 6.625 


Ac= 0.0001791 
AH = 0.0437 m. 
Da = 0.01639 m. 
V = 0.0185 m/s. 


=0.00000962 m. 
HAH) = 0.000707 m. 
Wh =0.000919 m. 
5K = 0.0005 m/s 


The uncertainty 6Re/Re = 0.0622 (6.22 %) or Re = 315.65 +19.63 


2. Colburn j Factor Uncertainty 


^ Re.Pr*'^ 
and 

J 

where. 


For 9 ”= 1300.36 W/ml 
V7’(73-n) = 3.97‘’C. 
Nua^ = 8.94 
/ = 0.0151 


6g” = 43.32 W/m". 
8Vr=0.1 °C. 
8 /Vm^ = 0.625 
8> = 0.0014 


The uncertainty Sj/j = 0.093 or 9.3 % 
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